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Gas Dynamics Equations: 
Computation 

< : Sh ° ck waves ' v ° rtlclty waves ' and entr ° Py waves are fundamental dlscontmmty waves 

in nature and arise in supersonic or transonic gas flow, or from a very sudden release 
(explosion) of chemical, nuclear, electrical, radiation, or mechanical energy in a limited 
space. Tracking these discontinuities and their interactions, especially when and where 
new waves arise and interact in the motion of gases, is one of the main motivations for 

m ■ 
m , 

numerical computation for the gas dynamics equations. 

The fundamental equations governing the dynamics of gases are the compressible 
Euler equations, consisting of conservation laws of mass, momentum, and energy: 



d t p + V • m = 0, d t m + V • (— — — ) + Vp = 0, d t {pE) + V • (m{E + p)) = 0, (1) 

P 

where V is the gradient with respect to the space variable x e R d , p is the density, 
v G R d is the gas velocity with pv = m the momentum vector, p is the scalar pres- 
sure, and E = ||v| 2 + e(r,p) is the total energy with e the internal energy, a given 
function of (p, p) defined through thermodynamical relations. The notation a ® b de- 
notes the tensor product of two vectors. The other two thermodynamic variables are 
the temperature 9 and the entropy S. If (p, S) are chosen as the independent vari- 
ables, then the constitutive relations (e,p, 8) = (e(p, S),p(p, S),9(p, S)) are governed 
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by OdS = de + pd(~). For a polytropic gas, p = Rp9, e = c v 6, 7 = 1 + ^-, and 

p = p(p, S) = K P ^e s ^, e = -^-p-r-'e 5 ^ = (2) 

7 — 1 7 — 1 

where R, c v , and k are positive constants, respectively. System ([[]) is complemented by 
the Clausius inequality: d t (pa(S)) + V • (ma(S)) > in the sense of distributions for 
any a(S) G C 1 , a'(S) > 0, to identify physical shocks. 

The Euler equations for an isentropic gas take the simpler form: 

m (R) 1x1 

dtp + V -m = 0, 9 t m + V-( ) + Vp = 0, (3) 

p 

where p(p) = /«oP 7 with constants 7 > 1 and k > 0. 

These systems fit into the general form of hyperbolic conservation laws: 

d t u + V • f (u) = 0, u G R m , x G R d , (4) 

where f : R m — >■ (R m ) d is a nonlinear mapping. Besides ([T]) and ([3]), most of partial dif- 
ferential equations arising from physical or engineering science can be also formulated 
into form (j2J) or its variants, for example, with additional source terms or equations 
modeling the effects of dissipation, relaxation, memory, damping, dispersion, magne- 
tization, etc. Hyperbolicity of system fll]) requires that, for all £ G S^ 1 , the matrix 
(£• Vf (u)) mxm have m real eigenvalues Aj(u, £), j — 1, 2, • • • , to, and be diagonalizable. 

The main difficulty in calculating fluid flows with discontinuities is that it is 
very hard to predict, even in the process of a flow calculation, when and where new 
discontinuities arise and interact. Moreover, tracking the discontinuities, especially their 
interactions, is numerically burdensome (see [H El H21 IT5]). 

One of the efficient numerical approaches is shock capturing algorithms. Modern 
numerical ideas of shock capturing for computational fluid dynamics can date back to 
1944 when von Neumann first proposed a new numerical method, a centered difference 
scheme, to treat the hydrodynamical shock problem, for which numerical calculations 
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showed oscillations on mesh scale (see Lax [H]). von Neumann's dream of capturing 
shocks was first realized when von Neumann and Richtmyer [23] in 1950 introduced the 
ingenious idea of adding a numerical viscous term of the same size as the truncation 
error into the hydrodynamic equations. Their numerical viscosity guarantees that the 
scheme is consistent with the Clausius inequality, i.e., the entropy inequality. The 
shock jump conditions, the Rankine-Hugoniot jump conditions, are satisfied, provided 
that the Euler equations of gas dynamics are discretized in conservation form. Then 
oscillations were eliminated by the judicious use of the artificial viscosity; solutions 
constructed by this method converge uniformly except in a neighborhood of shocks, 
where they remain bounded and are spread out over a few mesh intervals. 

Related analytical idea of shock capturing, vanishing viscosity methods, is quite 
old. For example, there are some hints about the idea of regarding inviscid g ases as 
viscous gases with vanishingly small viscosity in the seminal paper by Stokes (1848), as 
well as the important contributions of Rankine (1870), Hugoniot (1989), and Rayleigh 
(1910). See Dafermos [6] for the details. 

The main challenge in designing shock capturing numerical algorithms is that 
weak solutions are not unique; and the numerical schemes should be consistent with 
the Clausius inequality, the entropy inequality. Excellent numerical schemes should 
be also numerically simple, robust, fast, and low cost, and have sharp oscillation-free 
resolutions and high accuracy in domains where the solution is smooth. It is also 
desirable that the schemes capture vortex sheets, vorticity waves, and entropy waves, 
and are coordinate invariant, among others. 

For the one- dimensional case, examples of success include the Lax-Friedrichs 
scheme (1954), the Glimm scheme (1965), the Godunov scheme (1959) and related high 
order schemes; for example, van Leer's MUSCL (1981), Colella-Wooward's PPM (1984), 
Harten-Engquist-Osher-Chakravarthy's ENO (1987), the more recent WENO (1994, 
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1996), and the Lax-Wendroff scheme (1960) and its two-step version, the Richtmyer 
scheme (1967) and the MacCormick scheme (1969). See H El El HH CE1 HH EI] and 
the references cited therein. 

For the mult i- dimensional case, one direct approach is to generalize directly 
the one-dimensional methods to solve multi-dimensional problems; such an approach 
has led several useful numerical methods including semi-discrete methods and Strang's 
dimension- dimension splitting methods. 

Observe that mult i- dimensional effects do play a significant role in the behav- 
ior of the solution locally, and the approach that only solves one-dimensional Rie- 
mann problems in the coordinate directions clearly lacks the use of all the multi- 
dimensional information. The development of fully multi-dimensional methods requires 
a good mathematical theory to understand the multi-dimensional behavior of entropy 
solutions; current efforts in this direction include using more information about the 
multi-dimensional behavior of solutions, determining the direction of primary wave 
propagation and employing wave propagation in other directions, and using transport 
techniques, upwind techniques, finite volume techniques, relaxation techniques, and ki- 
netic techniques from the microscopic level. See [21 [131 El E0] • Also see [8l \TU \ [TT ] [T6" l l2Tj 
and the references cited therein. 

Other useful methods to calculate sharp fronts for gas dynamics equations in- 
clude front-tracking algorithms [51 E], level set methods [321 EE], among others. 
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